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^ : ABSTRACT 
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<^ | We determine the extinction curve in the zi = 0.83 lens galaxy of the grav- 

00 \ itational lens SBS0909+532 from the wavelength dependence of the flux ratio 

between the lensed quasar images (z s = 1.38) from 3400 to 9200A. It is the first 
measurement of an extinction curve at a cosmological distance of comparable 

^ ; quality to those obtained within the Galaxy. The extinction curve has a strong 

2175A feature, a noteworthy fact because it has been weak or non-existent in 

O ' most estimates of extinction curves outside the Galaxy. The extinction curve is 

fitted well by a standard Ry = 2.1 ± 0.9 Galactic extinction curve. If we assume 
standard Galactic extinction laws, the estimated dust redshift of z = 0.88 ± 0.02 
Oh' is in good agreement with the spectroscopic redshift of the lens galaxy. The 

widespread assumption that SMC extinction curves are more appropriate mod- 
1 els for cosmological dust may be incorrect. 
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1. INTRODUCTION 

Just as the properties of galaxies and stars evolve with cosmic time, so must the prop- 
erties of the interstellar medium. In particular, with changes in metallicity and ionizing 
backgrounds we might expect the properties of dust to change with redshift. Locally we can 
make precise measurements of extinction curves as a probe of the dust grains by comparing 
the spectral energy distributions (SED) of reddened and unreddened stars of the same spec- 
tral type (the pair method, Massa, Savage & Fitzpatrick (1983)), but the need to measure 
the spectra of individual stars limits this approach to the Galaxy, its satellites and its nearest 
neighbors. 

The average Galactic extinction curve is well modeled over a broad wavelength range as a 
function of only the ratio of total to selective extinction, R v (e.g. Cardelli, Clayton & Mathis 
1989, hereafter CCM; Fitzpatrick 1999). These extinction curves include a spectral feature, 
the 2175A bump, and appear to apply to a wide range of interstellar dust environments in 
the Galaxy. It is not, however, a universally valid law. While some lines of sight in the 
Magellanic clouds show similar extinction laws, others, like the 30 Doradus region of the 
Large Magellanic Cloud (LMC) show weaker 2175A features and steeper far UV extinction 
curves (e.g., Fitzpatrick 1986). The most extreme deviations from the Galactic extinction 
curves are found in the bar of the Small Magellanic Cloud (SMC, Gordon & Clayton 1998). 
Some newer studies (e.g., Clayton, Gordon & Wolff 2000) have found lines of sight with 
weak 2175A features and steep far UV extinctions in the Galaxy as well. M31 is the most 
distant galaxy to which the pair method can be applied; the M31 extinction curve seems 
similar to that of the Galaxy but may have a weaker 2175A bump (Bianchi et al. 1996). 

The studies of extinction in other galaxies are in many cases based on the examination 
of the attenuation induced by dust in the SED of the galaxies. But the physical properties of 
dust and its geometrical distribution affect in a complex way the intrinsic SED of the galaxy 
stellar population, which must be also modeled. The extinction curves inferred from these 
indirect methods are strongly model-dependent and, in many cases, not unique due to model 
degeneracies (Bianchi et al. 1996; Gordon, Calzetti & Witt 1997). In M33, Gordon et al. 
(1999) noticed that a depression found in the SED should be related to the 2175A feature, 
according to a radiative transfer model including attenuation by dust. A very weak 2175A 
feature is also found by Rosa & Benvenuti (1994) by comparing spectra of ionizing star 
clusters in HII regions of M101 with population synthesis models. However, the modeling of 
the SED of 30 starburst galaxies (Gordon, Calzetti & Witt 1997) implies that the dust has 
extinction properties similar to that of the SMC, lacking a 2175A feature. 

Some efforts have also been made in AGNs to study the effects of extinction by compar- 
ing the observed SEDs with conveniently selected templates. Crenshaw et al. (2001) have 
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compared the spectra of the Seyfert 1 galaxies NGC 3227 and NGC 4151 (assumed unred- 
dened), and found that the extinction in NGC 3227 increases very steeply toward the far 
UV (steeper than in the SMC) and that there is no evidence of the 2175A bump. In QSOs, 
the comparison of the observed spectra with composite spectra has not revealed traces of 
the 2175A feature (Pitman & Clayton 2000; Vanden Berk 2001). All the above-mentioned 
studies suggest that the presence of a significant 2175A bump in the average extinction curve 
of the Milky Way is rare (Pitman & Clayton 2000). In any case, the validity of this or any 
other conclusion about the behavior of the extinction beyond the Local Group is limited 
by the lack of a direct method to obtain the extinction curve. This limitation has serious 
consequences for the study of galaxies (see Gordon et al. 2000, and references therein) and 
in cosmology (see Falco, Kochanek & Muhoz 1998, and references therein). 

There are now two techniques which can extend the measurement of precise extinction 
curves into the realm of cosmology. One is based on the use of the relationship between 
the variations in color and absolute magnitude in SNIa light curves (Riess, Press & Kirsh- 
ner 1996). This technique has been used to estimate several reddening ratios from nearby 
supernovae which are in agreement with the MW average extinction curve. The other is 
an extension of the standard pair method to obtain extinction curves in gravitational lens 
galaxies. The method compares two different images of the same lensed QSO which follow 
different paths through the lens galaxy (Nadeau et al. 1991). As in the pair method each 
image has the same intrinsic spectrum but suffers from a different amount of extinction along 
each path through the lens galaxy. Using broad-band optical and NIR photometry, Falco 
et al. (1999) measured differential extinctions in 23 gravitational lens galaxies. In spite of 
the coarse spectral sampling, Falco et al. (1999) could fit the average Galactic extinction 
law to several lens galaxies, finding values for Ry in the range 0.64 to 7.2. This range in- 
cludes some extreme values for Ry whose physical significance is not obvious as the CCM 
law was inferred from values of Ry in the 2.5 to 5.5 range. Falco et al. (1999) also derived 
relatively accurate dust redshifts following the technique proposed by Jean & Surdej (1998). 
These studies relied on broad band photometry, which can lead to problems distinguishing 
between extinction and gravitational microlensing and cannot provide complete extinction 
curves over broad wavelength ranges. These problems can be overcome using spectra of the 
individual images, particularly using integral field spectroscopy (IFS) to obtain a 2D array 
of spectra of the lens and its images (Mediavilla et al. 1998). 

In this article, we use IFS to search for the 2175A feature in a cosmologically distant 
galaxy by selecting a gravitational lens where the lens galaxy shows significant extinction 
and the redshifted wavelength of the 2175A feature, 2175(1 + zi)A is below the Lyman break 
in the spectrum of the source quasar at 912(1 + z s )A. We performed the experiment on the 
two-image, l'.'ll separation quasar lens SBS 0909+532 (Kochanek et al. 1997), where the 
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z\ = 0.83 lens and z s = 1.38 source redshifts (Oscoz et al. 1997; Lubin et al. 2000) satisfy the 
wavelength criterion, and Falco et al. (1999) found there was significant differential extinction 
AE(B — V) = 0.20 ± 0.03 based on the image flux ratios measured in broad band imaging. 

2. DATA ANALYSIS 

2.1. Observations and Data Reduction 

The spectra of SBS 0909+532 were taken on January 18, 2001 with the 2D-fiber system 
INTEGRAL (Arribas et al. 1998), at the 4.2-m William Herschel Telescope (WHT). INTE- 
GRAL links the Nasmyth focus of the WHT with the WYFFOS spectrograph. We used the 
SB1 fiber bundle that consists of a rectangular array of 175 fibers, each 0'.'45 in diameter, 
covering an area of 7'/9 x &!7 on the sky, and a ring with diameter 90" formed by 30 fibers 
of the same size as those in the inner array. The distance between adjacent fibers is 0'.'5 in 
the central rectangle. At the entrance of the spectrograph, the fibers are aligned to form a 
pseudo-slit. 

We took three exposures of 1800s each. We used the 3001ines/mm grating, R300B, that 
yields a linear dispersion of 6. OA pixel -1 covering the wavelength range from 3400 to 9200A. 
The seeing, measured from a guide star was ~ 1". 

For our data reduction, we used the IRAF 7 software package. The procedure included 
bias subtraction, cosmic-ray rejection, spectra extraction, throughput correction, wavelength 
calibration, and sky subtraction. All these steps were applied to each one of the three 
individual frames and to an average frame obtained by combining them. For additional 
details on reduction and analysis of data taken with optical fibers, see Arribas, Mediavilla 
& Rasilla (1991). 

In Figure 1 we show the 2D distribution of spectra from the average frame over a spectral 
range from 6000 to 7200 A including the MgIIA2798 emission line. The positions of the two 
compact components A and B are also shown. 



7 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National 
Science Foundation 
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2.2. Continuum Maps 

From the 2D collection of spectra corresponding to the averaged frame (a subset of 
them around components A and B is shown in Figure 1) we interpolated 20 continuum maps 
covering all the observed spectral range. The central wavelengths and widths of the continua 
are included in Table 1. We fitted and removed the strong emission lines present in the QSO 
spectra. 

In Figure 2 we present six of these continuum maps corresponding to the averaged frame 
centered at 8727A, 7233A, 5739A, 4447A, 4006A, and 3487A. The two compact components 
of SBS 0909+532 appear well separated in the reddest map (8727A, seeing ~ 1'.'05). Notice 
the continuous fading of image B from the red to the blue maps. In fact, in the bluest map 
(3487A) the B image is no longer visible. This decrease in the B/A flux ratio towards the 
blue is explained by dust absorption in the lensing galaxy (Falco et al. 1999). 

2.3. A and B Spectra 

Because of differential atmospheric refraction (DAR), the location of any of the com- 
ponents in the telescope focal plane changes with wavelength. In fact, we measured a dis- 
placement of about 0'.'6 (greater than the diameter of a fiber face) between the bluest and 
reddest continuum maps. Although the displacement is modest, it implies that none of the 
spectra of Figure 1 can be exactly identified with one of the components throughout the 
full spectral range. This is a general problem of spectroscopy in a large spectral range that 
admits an a posteriori correction when 2D spectroscopy is available (Arribas et al. 1999). 
Following the procedure described in this reference, we calculated a DAR model from the 
location of component A in the 20 continuum maps distributed along the full spectral range. 
Next, from the original spectra we interpolated a new 2D set of spectra corrected for DAR 
(see Arribas et al. 1999 for the details). From these spectra we obtained a map in the 3500- 
9000A spectral range. We located the centroid of the A component in this map and used 
the coordinates offset of B from A given in Lehar et al. (2000) to infer the location of the B 
component. Finally, we interpolated at these locations the A and B spectra (see Figure 3) 
from the ensemble of spectra corrected for DAR. The changes from the uncorrected spectra 
are modest, as expected. 
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3. RESULTS 

3.1. Differential Extinction Curve from Continuum Images 

To derive the B/A flux ratios from the ensemble of continuum maps covering the full 
observed spectral range, we fit the maps using PSFs derived from the continuum images of 
a star observed on the same night. We used as fitting parameters the separation between 
the two QSO images, their relative intensity, the sky background, and the blurring of the 
PSF due to the seeing. We followed this procedure with each one the three independent 
exposures. After subtracting the quasar model, the pattern of residuals show no emission 
from the galaxy at any wavelength, so we conclude that our photometric model with only two 
point-like sources reproduces successfully the observed continua. Acording to the HST data 
(Lehar et al. 2000) the lens galaxy is not detected in the V band and its integrated magnitude 
is only 5% of component A flux in the I band. The location of the two quasar components 
determined by our photometric model is in excellent agreement with the results obtained 
with the HST by Lehar et al. (2000) which strengthens the validity of our photometric 
decomposition. 

In Table 1 and Figure 4 we present the mean magnitude differences, mg — tha, derived 
from the PSF fitting. To estimate the errors we have used the dispersion between the results 
obtained from each independent exposure. We include in this Figure the mg-mA differences 
from Lehar et al. (2000) and Kochanek et al. (1997) corresponding to broad band data; the 
agreement is excellent. 

3.2. Differential Extinction Curve from Emission Line Images 

Proceeding in the same way as for the continua, we have also computed the tub — rriA 
differences for the emission in the MgIIA2798 and CIII]A1909 lines. To estimate and subtract 
the continuum below each line, we fitted a straight line to adjacent points. The magnitude 
differences in the emission lines are approximately 0.3 mag above that of the continuum 
{tub — mA, see Figure 4). These differences are probably created by microlensing from the 
stars in the lens galaxy which can differentially magnify the continuum and emission line 
regions of the quasar because of their differing physical sizes (Schneider, Ehlers & Falco 
1988). While there are too few emission lines to estimate an extinction curve, the emission 
lines are consistent with the extinction properties estimated from the continuum. We have 
also tried to compute the m# — rriA difference from the CIV]A1548, but the maps were too 
noisy and the results of the PSF fitting were not consistent. 
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Table 1. 



A a 


W b 




A(m B -m A ) d 


3487 


252 


1.40 


0.20 


3565 


96 


1.45 


0.06 


3661 


96 


1.59 


0.05 


3757 


96 


1.67 


0.12 


3853 


96 


1.83 


0.10 


3949 


96 


1.77 


0.13 


4063 


96 


1.75 


0.10 


4159 


96 


1.63 


0.17 


4255 


96 


1.52 


0.08 


4351 


96 


1.52 


0.08 


4447 


96 


1.35 


0.09 


4743 


498 


1.07 


0.08 


5241 


498 


0.91 


0.03 


5739 


498 


0.83 


0.02 


6237 


498 


0.69 


0.02 


6735 


498 


0.55 


0.03 


7233 


498 


0.47 


0.03 


7731 


498 


0.43 


0.03 


8229 


498 


0.37 


0.02 


8727 


498 


0.37 


0.02 



a Contral wavelength of the continuum 
maps (A) 

b Wavelength bin of the continuum maps 

(A) 

c Averagc magnitude differences 
d Error in magnitude differences 
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3.3. Differential Extinction Curve from A and B Uncontaminated Spectra 

The method followed in Section 3.1 could be repeated with narrower continua to obtain 
the differential extinction curve with improved spectral sampling that could even match the 
spectral resolution of the data. However, the spectral sampling that can be achieved with 
this method is, in practice, limited by the noise in the data and the PSF. 

To improve the spectral resolution we can try, alternatively, and only because the lens 
galaxy is undetected, to obtain the extinction curve directly from the A and B spectra after 
correcting them for the cross-contamination induced by the seeing. We can estimate the 
contribution of A to the spectrum measured at the position of B by measuring the spectrum 
at position C located at the same distance from A as in the image B but in the opposite 
direction (see Figure 1). Similarly we obtain spectrum D as an estimate of the contribution 
to A from B. Our estimate of the seeing-corrected spectra are the differences A u = A — D 
and B u = B — C . 

We can check the results by comparing the spectral ratio to the estimates made by 
directly fitting the continuum maps. They agree reasonably well, with a ~ 0.15 mag drift 
between the blue and the red wavelengths. The drift appears to be due to small asymmetries 
in the PSF. We can improve the agreement by defining A u = A — (1 — e)D and B u = 
B — (1 + e)C with e = 0.1 to correct for the asymmetries. Figure 5 shows our final estimate of 
the magnitude differences [mg„ — toaJ(1/Ao) after smoothing by a a = 30A Gaussian. This 
curve is in good agreement with the difference magnitudes obtained from the continuum 
images at all wavelengths except in the vicinity of the emission lines. As stated in the 
previous section, this could arise from microlensing that is not affecting the lines to the 
same degree as the continuum. This implies that, in some spectral regions (like the one 
including the CIVA1548 line in SBS 0909+532) observations made with narrow-band filters 
could be problematic for obtaining the extinction curves since the contribution from the 
continuum and the strong emission lines must be separated. (Notice also the presence of the 
Mg IIAA2796,2803 absorption doublet arising from the lens galaxy). 



Assuming that any microlensing effects are achromatic, the magnitude difference (Fig- 
ure 4) is a simple function of the magnification ratio AM in magnitudes, the differential 
extinction AE(B — V) and the extinction law R(X) (Falco et al. 1999), 



4. 



DISCUSSION 




(1) 
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We fitted the m\ — rriB data from the continua using equation (1) and the CCM pa- 
rameterization of the extinction curve R(\). In a first step, the redshift is restricted to the 
value determined directly from the absorption lines of the lens galaxy (z = 0.83). The 
best fit is obtained for Ry = 2.1 ± 0.9, AM = —0.2 ± 0.2 and differential extinction 
AE(B — V) — 0.21 ± 0.02. A visual inspection of the fit (see Figure 4) shows reason- 
able agreement with the CCM extinction curve for the MW (x 2 /dof ~ 2). In spite of the 
uncertainties in R V) the fit of the 2175A bump is good in both shape and central wave- 
length. When, in a second step, changes not only in R v but also in z were allowed, we found 
that Ry was very poorly determined, but we obtained a useful estimate for the redshift, 
z = 0.88 ± 0.02. The close agreement with the spectroscopic redshift of the lens galaxy 
(z = 0.83, Oscoz et al. 1997; Lubin et al. 2000) demonstrates the possibilities of the dust- 
redshift technique (Jean & Surdej 1998; Falco et al. 1999) and reinforces the reliability of 
the extinction law. 

The detection of a significant 2175A feature in the extinction curve of a galaxy at 
z = 0.83 is very noteworthy. Many extragalactic extinction studies have suggested that the 
2175A feature is rare (see Pitman & Clayton 2000, and references therein). However, these 
studies were based mainly on starburst galaxies and AGNs, where the dust is subject to 
processing by radiation and shocks, similar to the environments of star formation regions of 
the LMC or the SMC, where the 2175A bump is weak. Thus, the detection of this feature in 
the extinction curve of the SBS 0909+532 lens (a normal early-type galaxy) and its absence 
in active galaxies follows naturally under the hypothesis of a connection between activity 
and suppression of the 2 175 A feature. 

It has been proposed that the presence of the 2175A feature is anti-correlated with 
the FUV extinction, with the feature weakening as the FUV extinction steepens (Clayton, 
Gordon & Wolff 2000). Accordingly, the rather low value Ry = 2.1 (although not rare in 
the extragalactic domain, see Falco et al. 1999 and references therein) that implies a steep 
FUV rise would be inconsistent with the strong bump (e.g. HD 210121, Clayton, Gordon 
& Wolff 2000). Within the uncertainty of the determination of Ry the average value for 
the MW (Ry = 3.1) would also be acceptable and in better accordance with the presence 
of a significant bump. Additional UV and IR data would be needed to better answer this 
question. 

5. CONCLUSIONS 

The application of the standard pair method to spectroscopic observations of the gravi- 
tationally lensed QSO SBS 0909+532 has allowed us to obtain the first determination of an 
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UV extinction law beyond the Local Group. Our conclusions are: 

1. We have detected a significant 2175A bump in a z = 0.83 galaxy. So far, most studies of 
extinction in the extragalactic domain have been based on starburst and AGN galaxies 
which seem to exhibit SMC-like dust properties lacking a 2175A feature. However, the 
lens galaxy of SBS 0909+532 seems to contain dust like that of the Galaxy. 

2. The Cardelli, Clayton & Mathis (1989) law provides a reasonable fit to the SBS 0909+532 
extinction curve with R v = 2.1 ± 0.9. This value includes, at ~ la, the average MW 
law (Ry = 3.1). However, values of Ry < 2.1 will be difficult to reconcile with the 
proposed experimental correlation between weak bumps and steep FUV extinction. 

3. The lens redshift estimated from the dust, z dust = 0.88 ± 0.02, agrees well with the 
redshift inferred from the absorption lines of the lens galaxy. 

We know of about 70 gravitational lenses, many of which show signs of differential 
extinction (Falco et al. 1999). By applying our techniques to a wider sample of these lenses 
we should be able to provide a broad outline of the properties and potentially the evolution 
of extinction curves with cosmic time. 
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Fig. 1. — The 2D distribution of spectra around the A and B components, in the 6000 - 
7200A wavelength range. The locations of the two QSO images are marked A and B. C and 
D mark the locations where two spectra were interpolated to correct the effects of seeing in 
the A and B spectra (see section 3.3) 
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Fig. 2. — The plots shows 6 of the 20 continuum maps described in the text. From top to 
bottom, left to right, the continuum central wavelengths are 8727A, 7233A, 5739A, 4447A, 
4006A, and 3487A. 



- 15 - 




Fig. 3. — Spectra corresponding to the A and B QSO images corrected for the detector 
response. The small box shows the original spectra. 
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Fig. 4. — Magnitude differences m# — uia obtained from continuum images as a function of 
inverse wavelength in the lens galaxy rest frame (filled squares). Open squares correspond to 
the B, V, R, I broad-band observations (Kochanek et al. 1997; Lehar et al. 2000). The solid 
line is the best fit of the analytical average MW extinction law to the observed extinction 
curve (filled squares) for a fixed z = 0.83. The values Ry = 2.1 ± 0.9 and differential 
extinction E(B — V) = 0.21 ± 0.02 were obtained in the fit. The filled triangles correspond 
to the m B — rriA differences computed from the emission in the Mgll and CIII] lines. The 
dashed line corresponds to the best fit for R v = 3.1. 
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Fig. 5. — The continuous line is the m Bu — rriA u magnitude difference curve obtained from 
the A and B spectra. The abscissa is the inverse wavelength at the lens galaxy rest frame. 
The squares joined by straight lines correspond to the data presented in Figure 4. 



